Universality classes of dense polymers and conformal sigma models 
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In the usual statistical model of a dense polymer (a single space-filling loop on a lattice) in two dimensions 
the loop does not cross itself. We modify this by including intersections in which three lines can cross at the 
same point, with some statistical weight w per crossing. We show that our model describes a line of critical 
theories with continuously-varying exponents depending on w, described by a conformally-invariant non-linear 
sigma model with varying coupling constant g 2 , > 0. For the boundary critical behavior, or the model defined 
in a strip, we propose an exact formula for the £-\eg exponents, he. = g 2 £{£ — 2)/8, which is shown numerically 
to hold very well. 

PACS numbers: 05.50.+q, 05.20.-y 
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Loop models are ubiquitous in low dimensional statistical 
mechanics, and have been studied for decades [1]. They have 
recently grown to play a major role in topological quantum 
computing J2]. 

Most loop models studied so far have to forbid intersections 
to be solvable. Their critical exponents can then be calculated 
using techniques of conformal field theoryj3t], Coulomb gas, 
or stochastic Loewner evolution (SLE) 14j,|5[]. 

Self avoiding walks, whose long distance properties de- 
scribe real polymers at interfaces, are the simplest of all loop 
models. In the ordinary, so-called dilute case, it is known that 
allowing intersections (hence obtaining self-avoiding trails) 
does not change the long-distance properties @]. The dense 
case, where a single self avoiding loop on a lattice is forced 
to occupy a finite fraction of the sites (and thus resembles a 
real polymer in a melt), is different. Allowing intersections 
does take the model to another universality class 0], which, 
however, shares many features with ordinary Brownian mo- 
tion, and does not seem to exhibit new families of critical ex- 
ponents. The first important result of this paper is that only 
allowing intersections where three lines cross simultaneously 
produces a very different behavior: a line of critical points is 
obtained, with central charge c = — 2, and continuously vary- 
ing critical exponents, 

Our second important result concerns the nature of this crit- 
ical line, which is very unusual in statistical mechanics. 

A convenient way to describe many loop models is to use 
a supersymmetric (SUSY) formulation, in which the degrees 
of freedom can take bosonic or fermionic values jjj, and the 
action is invariant under the action of a supergroup. The re- 
sulting field theories are however difficult to solve, in part due 
the lack of unitarity, and of current algebra symmetry. 

In the last few years, progress on one kind of such theories 
- cr-models on supergroups or supercosets - has been achieved 
in the framework of the AdS/CFT conjecture 

0E3I. An 

archetypal example is the principal chiral (PCM) model on 
PSL(2\2), which was found to be massless for a large range 
of values of the coupling constant g%. This is very different 



from what happens in ordinary groups, such as SU (2), where 
the PCM exhibits asymptotic freedom and spontaneous mass 
generation. The presence of conformal invariance and (su- 
pergroup symmetry is very interesting and potentially useful, 
yet, despite a lot of work, no complete solution of even the 
PSL(2\2) case, has been achieved full. 

We show in this paper that allowing crossings in dense 
polymers leads to close cousins of the PSL(n\n) models: 
cr-models on superprojective spaces (the super-analogs of or- 
dinary projective spaces) U(n\n)/U(n — l\n) x U(l). This 
identification has crucial consequences. It bridges the study of 
loop models with the one of a models, it gives direct access 
to properties of the cr-models both numerically, and, poten- 



tially, analytically using the techniques developed in 1112. fl3 1 . 
It also allows the determination of the critical exponents in the 
original geometrical problem. 

SUSY formalism for dense polymers. The universality 
class of dense polymers is generically obtained when one 
forces a finite number of self-avoiding loops or walks to fill up 
a fraction of space p > 0. When p = 1 for finite systems, one 
obtains Hamiltonian walks, the polymer limit of fully-packed 
loop models. For p = 1 the CFT is lattice dependent, i.e., 
universality breaks down Il4tl . Nevertheless, the loop model 
in Fig.Q]is in the generic dense polymer universality class. 

The CFT of dense polymers has c — —2; we review it be- 
low. It was discovered a few years back J3l that if one allows 
four-leg crossings the model flows to a different universality 
class with c = — 1, and trivial geometrical exponents. Such 
crossings imply that loops no longer conserve the lattice orien- 
tation of Fig.Q] indicating that a crucial symmetry is broken. 

To identify this symmetry we need to get into a bit of al- 
gebra. We consider the lattice in Fig. [T] and a transfer ma- 
trix T propagating vertically. We introduce a supersymmet- 
ric (SUSY) formulation lfl5i [HI]: each edge carries a Z 2 - 
graded vector space of dimensions m + n (resp. n) for the 
even bosonic (resp. odd fermionic) subspace (to + n, n > 
are integers). We label edges i = 0, 1, . . . , 2L— 1 for a system 
of width 2L. The Z 2 space is chosen as the fundamental □ of 
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FIG. 1 : Vertices, weights and sample configuration for dense poly- 
mers on a square lattice of width L = 3. Boundary conditions are 
free in the horizontal (space) direction and periodic in the vertical 
(imaginary time) direction. The alternating □, □ representations cor- 
respond to a lattice orientation, conserved along each loop. 



the Lie superalgebra gl(m + n | n) for i even (down arrow), and 
its dual □ for i odd (up arrow). T acts on the graded tensor 
product H = (□ <g> □) (this "Hilbert" space has in fact an 
indefinite inner product). 

To construct T for critical dense polymers we first observe 
that, for generic to, the tensor products □ ® □ and □ ® □ 
decompose as the direct sum of the singlet and the adjoint. 
The projectors on the singlet obey the Temperley-Lieb alge- 
bra relations Ef = mEi, [E.^Ej] = for \i - j\ > 2, 
and EiEi±\Ei = JSj (and here m = 0). The Ei can be 
expressed as quadratic terms in the SUSY generators, and 
are closely related with the Casimir 

d. We have T = 
TxT 3 • • • T 2L _ 3 T T 2 • • ■ T 2L _ 2 , where T, = 1 + xE,. By tak- 
ing either of the two terms in T, for each vertex, the expansion 
of Fig. Q] is obtained, with a power of x for each vertex, and 
a factor (n + m) — n = str 1 = m for each loop. The latter 
equals the supertrace in the fundamental representation (de- 
noted str) of 1, since states in H flow around the loop. This 
holds whether the loop be topologically nontrivial or homo- 
topic to a point. Isotropic dense polymers now correspond to 
in = and x = 1. Note that when m = 0, the tensor product 
□<£>□ is indecomposable; Ei can then be defined as the unique 
invariant coupling (on two sites) other than the identity. The 
rest of the discussion is unchanged. 

Letting x — > allows one to extract the spin chain hamilto- 
nian H oc — ^2 1 E% acting on TL\ the scale of H is chosen to 
ensure conformal invariance. The interaction is simply the in- 
variant quadratic coupling (Casimir), providing a natural gen- 
eralization of the Heisenberg chain to the gl(n + m\n) case. 

For such models, there is a corresponding continuum quan- 
tum field theory lllal . which is a nonlinear a-model with tar- 
get space the symmetry supergroup [here U(n + m\n)], mod- 
ulo the isotropy supergroup of the highest weight state (see 
|H EH for related non-SUSY examples, and E O for 
SUSY random fermion problems). Here we obtain U(n + 
m\n)/U{\) xU(n + m- l|n) = CP" +m - 1|n , a SUSY 
version of complex projective space. Moreover, the mapping 
shows that this model has a topological angle 8 = n. 

Dense polymers and sigma models. Let us now make 
things concrete: the fields can be represented by complex 



components z a (a = 1, . . . , n + to) and £ Q (a = 1, . . . , n), 
where z a is commuting, £" is anticommuting. In these coor- 
dinates, at each point in spacetime, the solutions to the con- 
straint z\z a + CaC" = 1 ( we use tne conjugation f that obeys 
(r)£y = for any rj, £), modulo U(l) phase transforma- 
tions z a i ► e lB z a , ( a i-> e lB C, a , parametrize cP" +m " 1|n . 
The Lagrangian density in 2D Euclidean spacetime is 

C = ^ [{d^-ia^zlid^ + ia^ (1) 

i6 

+ (dp - ia^didf, + ia^C] + 7T^ a " ~ d » a n)> 

Z7T 



where a, 
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+ CM a - {dz\)z a - (9Cl)C Q ] for 
fx = 1,2. The fields are subject to the constraint, and under 
the U(l) gauge invariance a M transforms as a gauge poten- 
tial; a gauge must be fixed in any calculation. This set-up 



is similar to the non-SUSY CP™" 1 model in BEO]. The 
coupling constants are <j£, the usual er-model coupling (there 
is only one such coupling, because the target supermanifold 
is a supersymmetric space, and hence the metric on the target 
space is unique up to a constant factor), and 6, the coefficient 
of the topological term (8 is defined modulo 27r). 

First we note a well-known important point about the SUSY 
models: the physics is the same for all n, in the following 
sense. For example, in the present model, correlation func- 
tions of operators that are local functions (possibly including 
derivatives) of components a < n\ + to, a < n\ for some n\ 
are equal for any n > m, due to cancellation of the "unused" 
even and odd index values. This can be seen in perturbation 
theory because the unused index values appear only in sum- 
mations over closed loops, and their contributions cancel, but 
is also true nonperturbatively (it can be shown in the lattice 
constructions we discuss below). In particular, the renormal- 
ization group (RG) flow of the coupling g 2 a is the same as for 
n = 0, a non-SUSY cr-model. For the case of CP" +m - 1|n , 
the perturbative /3-function is the same as for CP" 1 " 1 , namely 
(we will not be precise about the normalization of g 2 ) 



mg a 



(2) 



where £ = log L, with L the length scale at which the cou- 
pling is defined [see e.g. M22I1 . eq. (3.4)]. (The /3-function for 



is zero in perturbation theory, and that for g a is indepen- 
dent of 8.) For to > 0, if the coupling is weak at short length 
scales, then it flows to larger values at larger length scales. For 
8 it (mod 27r), the coupling becomes large, the U(n + m\ri) 
symmetry is restored, and the theory is massive. However, 
a transition is expected at 8 = it (mod 27r). For m > 2, 
this transition is believed to be first order, while it is second 
order for to < 2 lfl7ll . In the latter case, the system with 
8 = 7T flows to a conformally-invariant fixed-point theory. At 
the fixed point, a change in 8 is a relevant perturbation that 
makes the theory massive. 

For to — 0, the perturbative /3-function vanishes identi- 
cally. This can be seen either from direct calculations, which 



have been done to at least four-loop order 12211 . or from an ar- 
gument similar to that in |9|]: for n = 1, the cr-model reduces 



to the massless free fermion theory Il23ll C oc jfrC^iC 9/j,C an d 



further the 0-term becomes trivial in this case. Thus, for all cr- 
model couplings g 2 > 0, the n = 1 theory is non-interacting. 
The free-fermion theory is conformal with c = — 2, and 8 is a 
redundant perturbation, as it does not appear in the action (a 
similar argument appeared in Ref. 12m ). By the above argu- 
ment, conformal invariance with c = — 2 should hold for all n, 
and also for all g 2 and 6, though the action is no longer non- 
interacting in general. Thus the /3-function also vanishes non- 
perturbatively. In general, the scaling dimensions will vary 
with the coupling g 2 , so changing g 2 is an exactly marginal 
perturbation, though for n = 1 the coupling can be scaled 
away, so there is no dependence on the coupling in the expo- 
nents related to those multiplets of operators that survive at 
n = 1. Hence for n = 1, the exactly-marginal perturbation 
that changes g 2 is redundant. 

Introducing the six-leg crossings. For n — 1, the cr-model 
thus does not exhibit very interesting physics. It also de- 
scribes very few observables in the dense polymer problem. 
Indeed, the underlying algebra psl(l|l) does not admit any 
non-trivial invariant tensor, so the only £-leg operators present 
have £ = 0,2, and they are moreover degenerate — and part of 
an indecomposable block. These observables are expected to 
be present in all theories with n > 1 as well, and to not de- 
pend on the coupling constant g 2 . However, for n > 1, more 
observables are possible. E.g. £-\eg operators for all even i 
exist, and correspond to fully symmetric invariant tensors of 
psl(n|?i); there is no reason why the corresponding confor- 
mal dimensions should not depend on g 2 , and indeed we will 
shortly see that they do. 

For this, we need to be able to tune g 2 in the lattice model. 
We propose doing so by allowing not four-leg but six- leg 
crossings. This can be described most conveniently by go- 
ing to the hamiltonian formalism, and adding interactions that 
preserve the symmetry. Four-leg crossings would then trans- 
late into a perturbation of the type Pj,i+i which exchanges 
spaces at position i and i + 1. Since by construction our 
chain has alternating representations this is not possible within 
gl(n + m\n) symmetry, so forcing such crossings breaks the 
symmetry down to the orthosymplectic subgroup. On the 
other hand, six-leg crossings correspond to exchanging rep- 
resentations at position i,i + 2 while the one at i + 1 just goes 
through, and is perfectly compatible with the gl(n + m\n) 
symmetry (notice however that it breaks the extended sym- 
metry discussed in lfl2To . The hamiltonian then becomes 



Our first claim is that the continuum limit of (0) is de- 
scribed by the superprojective cr-model CP" -1 '" with cou- 
pling g 2 (w), at 8 = 7r. Note that we could more 
generally study the spectrum of the hamiltonian H oc 
- [ E % + wP i.i+2 + w 2 (E t E l+1 + E l+1 Ei)\. The sym- 
metries are unchanged, and one expects the continuum limit to 
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FIG. 2: Vertices and weights for dense polymers on the triangular lat- 
tice. When w = this is equivalent [25] to a Potts model with spins 
on the circles and arbitrary interactions within the gray triangles. 



be described by the same cr-model, with now g 2 (w,W2). This 
is confirmed by numerical calculations. Finally, a more pleas- 
ant realization of the same physics is provided by a model of 
dense polymers on the triangular lattice, where six-leg cross- 
ings can naturally take place; see Fig. |2] We shall call the 
Boltzmann weight of these vertices w as well, and the same 
conclusions will hold for this model as for the spin chain (f3j). 

We first check what happens for n = 1, where everything 
can be reformulated in terms of free fermion operators and 
their adjoints d /,,//, obeying {/«,/«,} = 0, {/<,/*} = 
(-1)*^ through Ei = (fj + f} +1 )(fi + f i+1 ) and P i>i+2 = 
(-iy + (//_! - - fi+i). Since both are quadratic 

it is easy to show that the continuum limit of <j3j is unchanged, 
with w only affecting the sound velocity and the fine structure 
of the Jordan blocks. 

One can easily argue that the ground state energy is the 
same for the n = 1 and n > 1 models, whence c = — 2 
independently of w. This is confirmed by transfer matrix cal- 
culations for the model in Fig. |2] The I = 2 exponent is 
conjugate to the fractal dimension of the loop, hence zero. 

Numerical study of the £ > 2 leg exponents then clearly 
shows that they are non-trivial, decreasing functions of w. 
To discuss this some more we place ourselves in the simplest 
case of free boundary conditions. The exponents at the special 
point w = are well known to be h® — /ii.i+f = e ( e ~ 2 ^ . We 
next assume that w — > oo corresponds to the weak-coupling 
limit of the cr-model, g 2 — * 0. This is qualitatively very 
reasonable: in the limit of large w, the system almost splits 
into two subsystems with gl(7i|n) symmetry involving only 
the fundamental or only its dual, with in both cases a simple 
interaction of the type P^j+i. Such models are well-known 
to be integrable, and their physics to be described by a weak- 
coupling limit not unlike the XXX ferromagnetic spin chain. 
In such a limit, we can analyze the spectrum using the minisu- 
perspace approach, that is, by analyzing quantum mechanics 
on the target manifold. The spectrum of the Laplacian on the 
ordinary projective space CP" 1-1 = U(m)/U(l)xU(m— 1) 
is well-known to be of the form Ei oc 41(1 + m — 1), so, set- 



(3) ting m = 0, we find that 12 



Here I in an 



integer, which we can identify using psl(n|n) representation 
theory with 1/2. Remarkably, /i™ c coincides with the known 
result h?f at w = (ordinary dense polymers) if we identify 
g 2 = 1 in that case. 

Conjecture for the exact exponents. We conjecture that 
the boundary conformal dimensions in our model are sim- 
ply linear in the Casimir of the associated representation of 



4 



-3.5 - 




FIG. 3: The two upper panels show the central charge as a function of 
the intersection weight w for width L strips of the triangular lattice. 
In the text we show analytically that c = —2 always. The lower left 
panel represents the effective coupling constant extracted from 
© using different values of I. The collapse on a single curve is 
quite striking. The lower right panel shows details of the exponents, 
in particular the region close to w = where convergence appears 
actually less good than on the previous curves. 



psl(n|?i). This is due to the structure of the perturbation theory 
where the vanishing of the dual Coxeter number — the Casimir 
in the adjoint — suggests exactness of the minisuperspace ap- 
proximation (see 12711 for a related case). A more thorough 
study of this perturbation theory, together with non pertur- 
bative arguments, will appear elsewhere 12811 . For now, we 
simply propose that the exponents be given by 



9o 



J(£-2) 



(4) 



where g\ is a decreasing function of w, equal to unity when 
w = 0, and vanishing at large w. 

This conjecture is compared with the results of exact diag- 
onalizations in Fig. [3] (lower panels), where we have repre- 
sented the function g^(w) as extracted from Q and various 
I. The different estimates collapse on a single curve over the 
whole range of w values, in agreement with the conjecture. 

For the model of Fig. [2] it is technically difficult to study 
operators with £ even. The er-model formalism can however 
be extended to £ odd, and the arguments leading to Q ex- 
tended to this case 112811 . Exact diagonalization of the spin 
chain hamiltonian on 2L = 18 sites yields results for £ even 
which look like the lower left panel of Fig. [3] except that w 
now has a different meaning. The sound velocity is deter- 
mined from analytical results for the n = 1 case. 

While the £-\eg exponents for the usual dense polymers 
(w = 0) agree with the general conjecture for g 2 a = 1, the 
fine structure of the spectrum at that point lfl3ll differs from 
the one of the sigma model. The situation seems similar to 
the one encountered in lL29ll for the supersphere sigma model, 
where the point w = is in fact singular. 

A related question concerns periodic boundary conditions. 
In this case, the known values of the bulk polymer exponents 



at w — are hg — = The fact that h e = 1 pro- 

vides an independent argument for the marginality of the w 
perturbation. Meanwhile, note that hi now do not have the 
minisuperspace form. For large w, one can however argue 
that the minisuperspace form remains valid, as is confirmed 
numerically. This suggests again that the point w = is sin- 
gular. It could also be that in the periodic case, the arguments 
that the minisuperspace should be exact for any w fail, which 
agrees with the expectations for a related model in 127]. More 
work is needed to clarify this point. 

We checked that, within numerical accuracy, staggering the 
chain produces similar results but with a coupling constant 
that now depends on w and the staggering parameter — that 
is, the 6 angle in the continuum limit. We e also studied the 
effect of coupling additional □ or □ representations on the 
boundary, which can be interpreted in terms of boundary 9 
angle 112811 . All the results are compatible with the er-model 
picture. 

In conclusion, we have shown that allowing intersections 
where three lines cross profoundly modifies the dense poly- 
mer problem. It gives rise to a critical line of conformal field 
theories, with central charge c = — 2, which can be identified 
with the long distance limit of a conformal sigma model such 
as those studied in the AdS/CFT correspondence . Our iden- 
tification leads moreover to the proposal of an exact formula 
(J4j> for the £-leg polymer exponents in the boundary case, and 
opens the way to tackling the sigma model using lattice tech- 
niques. 
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